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- SUMMARY 

An analysis of combined heat and mass transfer from a flat plate has 
been made in terms of Prandtl's simplified physical concept of the tur- 
bulent boundary layer. The results of the analysis show that for condi- 
tions of reasonably small heat and mass transfer, the ratio of the mass- 
and heat-transfer coefficients is dependent on the Reynolds number of the 
boundary layer, the Prandtl number of the medium of diffusion, and the 
Schmidt number of the diffusing fluid in the medium of diffusion. For 
the particular case of water evaporating into air, the ratio of mass- . 
transfer coefficient to heat-transfer coefficient is found to be slightly 
greater than unity. 


INTRODUCTION 

In recent years a considerable number of problems have arisen which 
involve the calculation of simultaneous mass and heat transfer at high 
speeds. The evaporative cooling of surfaces by air streams at high speed 
and at large Reynolds numbers is of considerable interest (ref. l). 

Although research on both mass and heat transfer has been conducted 
for many years and a number of analyses made, the status of the problem 
is such that it would seem desirable to present an analysis of the pro- 
cesses based on the simplified physical picture of the turbulent boundary 
layer and its laminar sublayer as originally conceived by Prandtl. 

The analysis presented herein is based on Donaldson's modification 
(ref. 2) of the Prandtl boundary-layer concept, which permits the calcula- 
tion of laminar -sub layer characteristics with a temperature variation 
through the boundary layer, and the use of Reynolds analogy in the tur- 
bulent region outside the sublayer. Such an analysis is, in general, 
subject to certain restrictive conditions since it must be assumed that 
the momentum boundary layer is unaffected by the mass- and heat-transfer 
processes. The application of the analysis is therefore limited to the 
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conditions of reasonably small temperature and partial pressure gradients 
across the boundary layer, small partial pressure of the material being 
transported in relation to the partial pressure of the medium of diffusion, 
and a Prandtl number of the medium of diffusion of order unity. Most 
gases, however, meet this requirement on Prandtl number. 

The requirements on the partial pressure of the fluids and the tem- 
perature gradient across the boundary layer are realized in at least 
several practical problems in aeronautics. In particular in the thermal 
de-icing of aircraft components, all the usual assumptions are met. 

The purpose of this paper is to present an analysis of the mass- and 
heat-transfer process in terms of a simplified physical picture of the 
turbulent boundary layer subject to the assumptions previously described. 


SYMBOLS 

The following symbols are used in this report: 

Cf skin-friction coefficient 

c specific heat of medium of diffusion at constant pressure, 

P Btu/ (lb) (°F) 

D diffusion coefficient, ft^/sec 

g acceleration due to gravity, ft/ sec^ 

k e mass-transfer coefficient or Stanton number for mass transfer 

kh heat-transfer coefficient or Stanton number for heat transfer 

m,k constants 

1 

n exponent of boundary- layer power-law profile, — = 

"0 

HCpg 

Pr Prandtl number, — jj — 

O 

Q heat transfer per unit area per unit time, Btu/ (sec) (fu ) 

Rl u l^l 

Re T Reynolds number, 

L Ml 

PqUqX 

Re v Reynolds number, 

X MO 
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Re S 


Reynolds number, 


Pp u Q 6 


r definitive ratio of total shear stress to laminar shear stress in 

boundary layer 

Sc Schmidt number, 

pD 

T absolute temperature, °R- 

u velocity, ft/sec 

W weight of medium diffusing from surface per unit area per unit 

time, lb/(ft^) (sec) 

x distance along surface from stagnation point, ft 

y distance normal to surface, ft 

5 boundary-layer thickness, ft 

laminar -sub layer thickness, ft 

concentration, weight of diffusing material per unit weight of 
medium of diffusion, Ib/lb 

K thermal conductivity of medium of diffusion, Btu/(ft ) (sec) (°F) 

P viscosity, slug/ (ft) (sec) 

V kinematic coefficient of viscosity, p/p, ft R /sec 

p mass density of medium of diffusion, slug/cu ft 

c mutual diffusivity, df^sion) 

(ft) (sec) 


.Subscripts : 

0 free -stream conditions 

adw adiabatic wall 

L conditions at edge- of laminar sublayer 


w 


wall 
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ANALYSIS 

Flow Over Flat Plate at Low Speeds 

The study of reference 2 relates the local skin-friction coefficient 
and the local Reynolds number based on boundary-layer thickness for the 
case of turbulent flow over a flat plate. The analysis described therein 
is based on the simplified physical concept of the boundary layer shown 
in figure 1. The boundary layer is assumed to be sharply divided into a 
turbulent region having a power-law velocity profile and a laminar region 
having a constant shear stress. The thickness of the laminar sublayer is 
given by the intersection of the turbulent -power -law velocity profile and 
the velocity profile of the laminar sublayer. This conception of the 
boundary layer is quite similar to that originally presented by Prandtl, 
which assumed a turbulent region and a laminar sublayer wherein the veloc- 
ity increased linearly with distance from the surface. 

The relations developed in reference 2 are in good agreement with the 
experimental evidence presented; hence, it may be assumed that for purposes 
of analysis the simplified physical picture of the boundary layer assumed 
by Donaldson (ref. 2) is justified. 

The general heat-transfer equation based on the temperature potential 
across the boundary layer, which is given in reference 3 is, in the nota- 
tion of this report, 

Q = k^gUQCp^ (T w - T 0 ) (1) 

Similarly, the mass-transfer equation based on the concentration 
potential across the boundary layer may be written as 


w = Vo^o (t W - V 


( 2 ) 


The heat transferred at the surface and throughout the laminar sub- 
layer results purely from conduction and may be written in the usual 
fashion as 


Q = - 


5 T 

K '<S7 


W 


St 

57 


( 3 ) 


Similarly, the mass transfer at the wall may be written as 


W = 


d 1 ! 

*w ^ 


w 



(*) 
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where o is the mutual diffusivity based on the diffusion characteris- 
tics of the substances involved in the process. In practical engineering 
terns, a is usually given as pDg where D is the diffusion coefficient 
for the particular process under consideration. 

The process of diffusion from the surface with its accompanying finite 
velocity of the escaping molecules is assumed not to affect the momentum 
boundary layer or heat transfer of the working medium. This assumption 
seems reasonable since the diffusing vapor represents only a small fraction 
of the medium of diffusion. 

The use of equation (4) (Fick*s law) in this form implies that the 
.partial pressure of the substance being transported is small compared with 
the pressure of the working medium. In addition, the mutual diffusivity 
o is defined for a field of uniform temperature; if large temperature 
differences exist; thermal diffusion will be superimposed over the mechani- 
cal diffusion and the simple form of Fick*s law will no longer apply. The 
mass -transfer analysis contained herein is therefore limited to cases where 
the temperature differences are small with respect to the absolute tempera- 
ture of the medium of diffusion and where the partial pressures of the 
diffusing material are small compared with partial pressure of the working 
medium. 


It is assumed that the temperature and concentration gradients at the 
edge of the laminar sublayer are given by 


ST ~ T L ~ ?w 
l 5 l 


( 5 ) 


St) ^ ~ ^w 

~^L " 5 L 


( 6 ) 


This assumption necessarily implies that the thickness of the laminar sub- 
layer is identical for both heat and mass transfer. There is, however, no 
a priori reason to expect that the physical thickness of the sublayer would 
•differ for the two processes. Donaldson (ref. 2) has shown that for the 
incompressible case with zero heat transfer, the laminar -sub layer thickness 
is dependent on a characteristic Reynolds number based on the velocity at 
the edge of the sublayer, the laminar-sublayer thickness, and the viscosity. 
Such a characteristic Reynolds number may be assumed applicable to heat and 
mass transfer provided the transfer rates are low. From equations (l), (3), 
and (5), and (2), (4), and (6), the following relations may be obtained: 


khCOgUo^O < T v " T 0) 



( 7 ) 
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, • _ , \ T W “ 

kePO^ ~ “ °L 8^ 


( 8 ) 


Rearranging terms and multiplying the left-hand side of equation (7) 
hy Mq/mq> K q/* 0-» and 6 / 6 yield the following 


PQ U 0 5 c p,QSHp K 0 V = T w - T L 


MO 


'0 "L 


•t T 5 T w - Tq 


hut 


c p, 06 P 0 = 

c o 


If— = Pr 0 


and 


p 0 U 0 5 p 

“o 8 


__ K 0 8 L T W " t L 
k h Re 6 Pro 6 “ t w - T Q 


(9) 


A similar result may he obtained for equation (8) hy rearranging 
term*, multiplying the left-hand side hy PqA’-O' 5 /h, and n L /n L and sub “ 

stituting = pjO^g. 


hut 


therefore 


M L MO = ^ 

^o pl d l 6 ~ \ - ^0 


MO Bl 

k e Reg Sc L — — 


% - 

Hw “ “HO 


( 10 ) 
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It is interesting to note the similarity between equations (9) and (10) 
with the exception that in equation (9) the Prandtl number is based on 
stream values and the Schmidt number in equation (10) is based on values 
at the edge of the laminar sublayer. In general, however, this is not 
important since the Prandtl number for air varies only slightly with tem- 
perature and the Schmidt number for many combinations of fluids is also 
only a slight function of temperature. 

rp — T ^ 

The temperature ratio — and the concentration ratio — 

T w “ T 0 \ 

may be evaluated in terms of the laminar-sublayer velocity u^ and the 
stream velocity u Q with the aid of Reynolds analogy. The following 
derivation is similar to that found in many heat-transfer text books } see, 
for example, reference 4. It is repeated here to obtain specific rela- 
tions required by the present analysis. 

Consider an area in the turbulent region of the boundary layer as 
shown in figure 2. Let 3 be a small mass of fluid which is transported 
per unit time, per unit area across the plane A-A. Suppose that 3 pene- 
trates upward through A-A to a region of higher velocity u', lower tem- 
perature T 1 , and lower concentration t] ' . In the steady-state condition 
(completely developed turbulence), an equal mass 3 must be transported 
downward across plane A-A to a region of lower velocity u, higher tem- 
perature T, and higher concentration t). 

The exchange in momentum across A-A is 3u - 3u'. This change in 
momentum must be balanced by a change in the shear stress such that the 
"virtual turbulent shear stress" must be equal to 3u - 3u*. 


T t = "P( u “ u' ) (11) 

Similarly the heat transported across A-A must be given by 

Qt = + 3 gc p (T - T') (12) 

and the mass transported across A-A by 

W t = + Pg( T l - tj») (13) 

To eliminate 3, which is unknown, from the equations, equations (12) 
and (13) are divided by equation (ll) as follows: 
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gc p (T - T') 


u - u' 


W + 


n-jT 


u - u 


r § 


From equations (3) to (6), it is found that 

?w - T L 


Q = 


‘V 


W - PL D L g 8- 




(14) 

(15) 


(16) 


(17) 


Also the shear stress at the edge of the laminar sublayer is given in 
reference 2 as 


T 


L " ^ &l 


(18) 


In accordance with the concept of a boundary layer with a sharp division 
between the laminar and turbulent regions, it is therefore assumed that 
heat and mass transfer in the turbulent region result entirely from tur- 
bulent interchange, that conductive effects are negligible, and further 
that in the laminar sublayer heat and mass transfer occur purely by con- 
duction and diffusion, respectively. 

Therefore at the edge of the laminar sublayer, = Q, W t = W, and 
T = and the following results. may be obtained from equations (14) 

. . T _ T . _ t l ' T 0. 

to (18) since 

K I 

L 5 l = _ s c p,l( t l - t q) 

U L U L " U 0 

" L5 l 


T w - T l SCp^L U L 


t l - T 0 


= Pr 


U L 

L u Q - u L 


u 0 - U L 
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Now 


plPl 


V - 





= _^L_ U L = Sc U L 
‘ ^0 P L°L u 0 " U L CL "0 “ U L 


X - T L 

T w - T l T w - T l T l - Tq 

T w - T 0 T L " T 0 + T w - T L T v " T L 

+ T L “ T 0 


U T 


U T 


Pr- 


L u 0 - U L 


Pr T — 
L u o 


Ut 


1 + Pr- 


■ L - 1 - — (1 - Pr L ) 

L ^0 " U L u 0 L 


Similarly, 


% - 

X " ^0 


U L 

u 0 


U T 


Ur 


(1 - Sc L ) 


Donaldson (ref. 2) has established that 


u 0 


n(r - 1) V L 


u 0 5 


1 

n+1 



n(r - 1) 


k 


2 


n 

ln+1 



(19) 


( 20 ) 


( 21 ) 


( 22 ) 
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where 


n(r - 1) 


k 


2 


is found to represent a characteristic Reynolds number 


of the laminar sublayer Re^ which has a constant value for a given value 
of n. A value of Re^ = 158 is found in reference 2 by a comparison of 
the analytical expression for the skin friction with an experimental drag 
law. 


C-p = 0.045 Reg 


-1/4 


This value of Re^ is in good agreement with the critical Reynolds number 
of the laminar sublayer as determined by Donaldson (ref. 2) from the re- 
sults of von Karman (ref. 5). Equations (21) and (22) may be rewritten 
as 


1 1 

u 0 


s 

now 

_ ^i/pl _ Pl Po 
y o " >V p o "Ho Pl 

but p-r / (Iq can be expressed to a close approximation as 

Ml _ /^\ m 
MO \ T 0/ 

Since the static pressure is constant throughout the boundary layer, 

PD _ 

. Pl T 0 
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Therefore, 


and 



"0 




(23) 


n n(l+m) 



(24) 


Solving equations (9) and (10) for k h and k e , respectively, and using 

the relations given by equations (19), (20), (23), and (24) give the 
following results: 


2 



2 



It is interesting to note that for the usual values of n=7, m=0.76 
and Re^ = 158 equation (25) becomes 
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Pr 

Pr 


- Re* -1 / 4 (0.0225) 


kh = 


^0 

*L 


1.32 




1 - 1,883 R6g 


T \°* 22 

%) (1 - Pr ^ 


For the case of low heat transfer, that is, low temperature differential 
across the boundary layer, 



Then 


k h = 


-1/4 

0.0225 Reg 

3178 


1 - 1.883 Reg 


(1 - Pr 0 ) 


For a flat plate, empirical results give 


Reg =0.37 (Re x ) 4 / 5 


Therefore 


0.0289 (Re x ) -1 / 5 

k h - -1/10 T \ 

1 - 2.13 Re x (1 - Pr 0 ) 

This equation is quite similar to the equation given on page 117 of refer- 
ence 4 and for air differs only slightly in the . value of ■ the constants in 
the numerator and denominator. Values of kh calculated from this equa- 
tion and from the equation given in reference 4 agree within several per- 
cent, the latter equation giving slightly lower values. This deviation 
is well within the experimental error found in most heat -transfer measure- 
ments. It can be seen, therefore, that equation (25) reduces essentially 
to a previously established law for heat transfer on a flat plate for the 
case where the temperature differential across the boundary layer is small 

The ratio of k e to kh as given by equations (25) and (26) is as 


follows: 
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Pro K 0 ^0 

In. the previous equation the term = — is obviously equal 

c Pr L k L 

to and hence the equation may be rewritten as follows: 

C P,L 



For the technically interesting case of water evaporating into air, 
the property values of both air and water are well known. In addition, 
the diffusion coefficient D for water vapor diffusing through air is the 
most reliable yet established. 

The limitation of the mass-transfer relation mentioned earlier, that 
is, low partial pressure of the diffusing material with respect to the 
pressure of the medium of diffusion, limits the applicability of equa- 
tion (27) to cases where the wall tempera-cures are much less than the 
boiling temperature of water. 

In addition, the analysis is valid only for the condition where the 
ratio Tq/t w is nearly unity and, hence, T^/Tq is also nearly unity. 

Since in equation (27) Tj/Tq is raised to a power considerably less than 
one and since the ratio c p,o/ c p,L is only slightly affected by tempera- 
ture variation, equation (27) may be written to a close approximation as 
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1 - 


k h 


1 - 


Rgj 

Rec 


1_ 

n+1 


(1 - Pr L ) 


Re L V 

Re S/ 


n+1 


(1 - 


Sc L ) 


(28) 


In the range of temperatures of practical interest (32° to 100° F), 
(for air) the Prandtl number is approximately equal to 0.71 and for water 
diffusing into air the Schmidt number is approximately 0.60. The varia- 
tion of with Reg is shown in figure 3 for the value of 

Re^ = 158 determined in reference 2 and several values of n. Also shown 

on the figure are curves calculated at a constant value of n equal to 7 
and values of Re-g equal to 100 and 300. The range of Reynolds numbers 

Reg shown correspond to a range of Reynolds numbers based on distance 

from the leading edge of 3.4X10^ to 1.91X10^. 

It is apparent from figure 3 that the ratio of the mass-transfer 
coefficient to the heat-transfer coefficient decreases slightly with in- 
creasing Reynolds number. The ratio also increases with increasing n, 
but both the effect of Reynolds number and n are quite small and it 
would be difficult to isolate either experimentally. The curves shown 
for values of Re^ equal to 100, 158, and 300 at a constant value of . 

n = 7 indicate that the effect of changes in Re^ is so small as to be 
negligible . 

The independence of the ratio k e /k^ with critical Reynolds number 
Re^ shows that the assumption made previously about the constancy of Re^ 
with or without heat transfer is not critical for purposes of the analysis. 

From an over-all standpoint, it would therefore be expected that the 
ratio kg/k^ can be taken as approximately 1.05 regardless of whether or 

not the flow is transitional with values of n of the order of 4 or fully 
turbulent with values of n equal to or greater than 7. 

Although there are few actual experimental data on the evaporation of 
water into an air stream from a flat surface, it has been found that the 
surface temperature of bodies wetted by a water film, calculated by 
assuming k e equal to k^, agrees quite well with experimental results 
(see, for example, ref. l). An unpublished experiment by Coles and Ruggeri 
on the mass and heat transfer from an iced flat surface for a wide range of 
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altitudes and Mach numbers shows that the ratio of the mass -transfer coeffi- 
cient to the heat-transfer coefficient is approximately 0.94. These results 
were obtained with rough, iced surfaces; consequently only qualitative agree- 
ment with the present analysis should be expected. 

It is evident that the analysis used herein lends itself to compara- 
tive results of heat and mass transfer but can be used to find absolute 
values of k e and k^ only when the value of n is known. For ex- 
tremely high Reynolds numbers, the value of n may be 10 or greater and 
this would greatly affect the values obtained from equations (25) and (26). 

The calculation of the ratio kg/k^ for fluids other than water into 
air can be accomplished easily provided a reliable value of diffusion 
coefficient is available. It is possible that certain fluids may exhibit 
properties such that the value of Tj^Tq could be significant without 

invalidating the assumptions of the analysis . In such a case, the general 
form of equation (27) must be used and Tj^/Tq obtained from a knowledge 
of T w /T 0 and trial -and-error solutions of equations (19) and (23). 


Extension to High-Speed Flow 

For high-speed flows wherein the frictional temperature rise in the 
boundary layer is appreciable, it is usual to write the heat -transfer 
equation in the following manner 


Q = k h p 0& a 0 c p (T w - T adw ) 

The use of the adiabatic wall temperature in place of the stream 
static temperature has been found to result in a satisfactory correlation 
of k h with Reynolds number for both low- and high-speed flows (up to 
Mach numbers of 2). 

The relation given by equation (25) would therefore be expected to 
hold over the whole speed range of current interest for aircraft icing 
given in reference 1 (Mach numbers from 0 to 1.5). The relation of equa- 
tion (26) is unaffected by speed provided the conditions stated previously 
with regard to the vapor pressures and vapor-pressure gradients in the 
boundary layer are still fulfilled. For the usual values of .interest in 
icing, these values are well within the acceptable limits. The ratio of 
k e /kh as given by equation (28) will therefore be independent of Mach 

number for the whole range of flight speeds of current interest. 
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CONCLUSIONS 

An analysis of combined heat and mass transfer from a flat plate has 
been made in terms of Prandtl's simplified physical concept of the bound- 
layer . The results of the analysis show that for conditions of r e a ~ 
sonably small heat and mass transfer, the ratio of the mass- and heat- 
transfer coefficients is dependent on the Reynolds number of the boundary 
layer the Prandtl number of the working fluid, and the Schmidt number of 
the diffusing material in the medium of diffusion. For the particular case 
of water evaporating into air, the ratio of mass-transfer coefficient to 
heat -transfer' coefficient is found to be slightly greater than unity. For 
the particular case of aircraft icing, it is shown that the results of the 
analysis are valid up to the maximum Mach number at which icing might 
occur, that is, 1.5. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 12, 1953 
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